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1. Introduction 

One of the important fields of modern ionospheric investigations is the study of 
small-scale irregularities in E- and F- layers of the ionosphere, that affect radiowave 
propagation and functionality of different HF and UHF radiotools. One of the most 
investigated types of irregularities is E- and F- layer irregularities produced as a result 
of growth of two-stream and gradient-drift instabilities. The theory of such instabilities 
is under development for a long time but still is not finished HI [191 ITT] . 

The usual condition for the growth of such irregularities is the requirement of 
different velocities of electrons and ions, most significant at altitudes 80-120km. At 
these heights ions are 'unmagnetized' - their motion is controlled by neutral component 
motion. At the same time the electrons motion is controlled by auxiliary electric and 
magnetic fields - electrons are 'magnetized' [ISl HI [19] . 

But it is clear that this requirement significantly limits the validity region of the 
current instabilities theories, that is why it is important to obtain a theory of ionacoustic 
instabilities without this limitation (see for example [H]). 



2. Basic equations 



2.1. Two- fluid magnetohydrodynamic equations 



As basic equations for obtaining dispersion relation we will use two-fluid magnetohy- 
drodynamics (MHD) equations in form [15] : 



where a = e,i, normalized elastic collision frequencies: 



(1) 



(2) 



Han is effective mass of charged particles during elastic collisions with neutrals (for 
electrons /ign is very close to electron mass, for ions /i^ could vary around half of ion 
mass) : 



rria + mr, 



(3) 



and z/^„ - elastic collision frequency of the ions (electrons) with neutrals. We also suppose 
here that the charged particles do not interact with each other through the collisions 
and interact only through electromagnetic field. We also exclude all the viscidity effects, 
that usually are not taken into account [T3]. In this case we take into consideration 
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(following to the [15]) only elastic collisions. In detail the approximations used are 
listed in Appendix A. Most of these approximations are valid at heights below 200km 
both for quiet and disturbed ionospheric conditions. 

It must also be noted, that we throw out a lot of terms from the MHD equations ([T]): 
ambient magnetic and electric fields are supposed to be constant; any magnetic 
field variations are not taken into consideration; gravitation field is not taken into 
consideration; recombination and ionization processes are not taken into consideration; 
neutral component motion is neglected. This allows us to neglect a lot of instabilities 
and effects (see, for example [1]) and simplify the analysis. 



2.2. Zero order solution - quasihomogeneous and static case 

Zero order approximation connects nondisturbed (quasihomogeneous and static) values 
of the particles density N^q, average motion speed V^aO; ambient electrical and 
magnetic fields and collision frequencies. When density, fields and collision 
frequencies are given, the zero order approximation defines average motion speed of 
charged particles in any point of space and time. 

As one can see from ([T]), the zero-order approximation is defined by the system: 

' V (iVaol^ao) = 

-m„Ar^o( V^ao V) V^«o - V(^aoA^ao)+ (4) 
= e(Z,Ar,o + AT^o) 



In the simplest case of weak velocity gradients, when we could neglect Lagrange 
term (V^aoV)'^oO; the system (H has a well known solution [TTf [15]: 



aO 



D 



Ta- 



where operators of diffusion Da, thermodiffusion Dxa and conductivity are: 
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(6) 

(7) 

(8) 
(9) 
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K = —X (10) 

ma 

D^a = —A^a (11) 

Oa = -Aa 12 

+ -I 

J. _ ^^(^) 

A. (13) 

and functions K„[x\Ki^[x) are tabulated (for example in |T7]) for taking into account 
not only MHD effects, but kinetic effects too. 

It is important to note that for our next consideration the exact expression (|5]) 
for zero-order solution is not very significant for us, and below we only suggest that 
the solution V^q,q = I/^q [Na-, Ta, ^o, ^o, ^ari) exists and is unambiguously determined 
by its arguments. So by the zero-order solution we mean an equation Q that defines 
an average motion speed \^ao as a function of ambient conditions and which could be 
solved analytically ([5 13) in simple cases or numerically in more complex cases. 



2.3. First order solution - nonstatic inhomogeneous case 

One of standard approaches to the MHD equations analysis is a geometrooptical (GO) 
approximation, the validity of which is defined by smallness of the parameter 

^^(^(^)) ^^^^ 

where k is irregularities wave vector and typical range of changes of parameter P 
(for example electron density). 

When the GO approximation is valid, the solution for small variations of the 
parameters N,\^ can be found in form: 

5Na{T',t) = e-*'^(^^'^*)/'^A/'„i(;ur,/it) 

sVai^, t) = e-^'^(^^''^*)/^ V^«i(/ir, jjt) (15) 
6^{-r',t) = -V (e-^^('^^'^*)/^$i(/i r',/it)) 

Geometrooptical phase ip^fi'f', fit)/ fi (or eikonal) for plane waves is related to wave 
vector k and complex frequency of the wave u + i'j hj the following definitions: 
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k 

LU + i'J 



- \/ij{fi't,fit)/^ 

di/j{fi'f', fit) 
fidt 



(16) 



The first approximation gave us the system of equations: 

Z,iV,,i = -iV,i-l$i(v^)2 
where, by taking into account the zero-order approximation (|4]): 



(17) 



Pal 
Pa2 
Pa3 



V^aOV^aO , .Jf ^ . .dtp 
hzKaOV^ + «^ 



Ta - UV^ + 



aO 



N, 



aO 



PaA = iZaeNao (vV') 

PaO^al = Pa5^al + al X Pa6 + Palpal 



Pa5 
-P 



-maN^O ( ^ ( ^ ) + ^ C^^aO V^) + ^an 



a6 



-ZaCNaoBo 
Palpal = -rria ( (A^ao'^^al V ) "^ao) 



(19) 

(20) 

(21) 
(22) 

(23) 

(24) 
(25) 



To make the following analysis easier, the system (17) is written in operator form, 



where operators P„ are matrix operators in partial derivatives over the eikonal if). It is 
clear that in this form the system looks pretty simple and solvable. 

3. Dispersion relation 

3.1. Obtaining the dispersion relation 



From (17) one can see that the system is linear and, in case of existence and uniqueness 
of the inverse operator P^^ (22) it can be solved. After excluding ai from (17) 



the equation connecting the density N^i and electric potential $i variations has the 
following form: 







(26) 
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(27) 



where coefficients are: 

and A{...) - an arbitrary function of arbitrary parameters that does not have zero values 
at the investigated region. 

Now we can recall that our plasma has two types of particles and its characteristics 
are defined by the system of equations: 

CH(^)iV,i + C2i(^)$l = 

CleWiVel + C2eW$l = (28) 

Z,iVa = -iVei-^$i(vV^)2 
Sometimes, for example, when analyzing thermal variations of electron density 



(that cause incoherent scattering), the self-coordinated term ^^i{\/ipy in (28) can 
not be neglected - scatterers size has order of Debye length and this term becomes 
significant. But in this very task we can neglect this term, following to many authors 
(see for example [TT]). 



It is clear that existence of solution of (28) is determined by consistency of these 



equations. The consistency condition in our case has the following form: 

CuWC2eW + Z,CM)C2^{'ip) = (29) 

It connects different partial derivatives over the eikonal ip with each other and can 
be referred as dispersion relation. As one can see, the dispersion relation has symmetrical 
(as it was expected earlier) form. 

3. 2. lAQV approximation 



It is clear that existence of dispersion relation and its exact form (29) depend on 
existence and properties of inverse operator P^o- As it has been shown in Appendix 
B, the inverse operator can be easily found in case when Lagrange term Py in 
Pao can be neglected. Below we call this approximation as 'Irregularities under 
approximation of quasihomogeneous velocity' (lAQV). As preliminary analysis has 
shown, this approximation is valid for wavenumbers 0.1-10 under most ionospheric 
conditions at altitudes below 200km and for variations of average parameters not faster 
than 100m (for faster changes the GO approximation becomes incorrect). 

In the lAQV approximation the inverse operator Pao has the following simple form: 



,y _ b (&7) Pi, + PlJ + P.,Pa,{h X 7) 

J T-> f T-,2 , T-,2 \ 



^aO J - P I p2 ^1 

Q5 l-T o6 "T aSj 

where b - unity vector in direction of 'Pae (and antiparallel to the magnetic field). 

It must be noted that lAQV approximation does not mean neglecting the Lagrange 
term (^T^av) '^a in basic equations (1), but only neglecting (^T^aiv) ao term in 



operator P^o (22), all the other terms are the same. 
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Summarizing, the dispersion relation in lAQV approximation has the following form 



(29, 18 24, 30) 



3. 3. The basic structure of the dispersion relation 

Let us briefly analyze the structure of the dispersion relation by defining function that 
does not have zeroes: 



A(...) = P«5 {Pl, + Pl) 



(31) 



This leads to the following coefficients of the dispersion relation (29): 



Pi,) + 



ClaW=^PalPa5 {P^e 

— (^a2bj (b^as) Pa6 ~ -Po5^a2^a3 + 
~PabPaQP a2{b X ^as) 

= - (bp 0.2) (bPaA) P^e - PlJ^a2'PaA + 

— Pa^PaQP a2{b X P 



(32) 



When taking into account ([T8|24) it becomes clear that coefficients (32) are 
polynomials over the and have the form: 



at J 



+ 



+-Rla2 (^) + -Rial (^) + PlaO 
= R2a2 (f )' + R2al (f ) + R 



(33) 



From this consideration it becomes clear that dispersion relation (29) is a 6th order 
polynomial over the and has no more than 6 solutions. 

3.4- Simple representation of coefficients 

To simplify the solution technique in homogeneous case, in the work p!3] a new complex 
variable was defined: 

dip 



dt 



+ 



In our inhomogeneous case we define the following new variables: 

1 



dip 

In 

N = V"^ + 



N 



t,fj, 
''anN 



N y aO 



In this case the operators (18] 24) become: 

Pal = i {<^aN) 



(34) 

(35) 

(36) 

(37) 
(38) 

(39) 
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P^2 = iN^ot*^ (40) 

pZ = -iTa^N (41) 
k -\- k ^ 



Pa4 = iZaCN, 



(42) 



Pa5 



-im, 



a6 — —ZaCNaoB Q 
Pae = ZaCNaoBo 



(43) 
(44) 
(45) 



where * is a complex conjugation. 



Summarizing, the dispersion relation has the form (29), where its coefficients are 
defined by ( p2p5p 5| 



4. Dispersion relation for weak gradients case 
4.I. Weak gradients approximation 

The vectors Pa2, Pas, Pai in dispersion relation are parallel in case of homogeneous 
ionosphere and not parallel in case of inhomogeneous ionosphere. After substituting 
Pa2-, Pa3, Pa4 (4o|42) iuto (32) and taking into account the properties of vector product 
we obtain the following: 



ClaW=PalPa5 {P^e + Pa,) + 



k iqh 



p2 I 



k 



N 



p2 



+ 



+i2NaoTaPa5Pa6l'fn 

C2.(^) = Z^eNl^ (Pl^ (hRe{tN)f) + 

+Z^eNl, [pl, (Re(tr,)f) + 

iZ^eNl^ [{Re{tN)) {lm{t m)) PJs) + 

' (bx ReCk 



(46) 



N 



+ 



N, 



p2 



The first two terms in each relation (46) are the terms that correspond both to the 



inhomogeneous and homogeneous dispersion relations, the last terms correspond only to 
changes of dispersion relation due to inhomogeneities presence. It is clear that in both 
cases (gradients are parallel to the magnetic field and gradients are perpendicular to the 
magnetic field) the difference between dispersion relation for homogeneous case and for 
inhomogeneous one does exist. But if the gradients are weak enough (or wavenumbers 
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are high enough) we can neglect the changes of dispersion relation and solve only 
simplified one, that corresponds to the homogeneous case: 



k Nb 



p2 I 



k 



N 



p2 



N, 



+ 



C2aW = Z^eNl^ (Pl^ (bRe{ k 
+Pl,[Re(t ^)f^ 

As one can see, this approximation is valid, when: 



(47) 
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» 1 



(4J 



(49) 



(50) 



The condition (48) is equivalent to the GO validity condition (14): 



N 



» 



The condition (49) is valid when analyzing scattering almost perpendicular to the 
magnetic field or when gradients are sufficiently small: 



k 



N 



a5 



>> 



(bRe{t 



N, 



a6 1 



(^Re(y 



(bt 



N 



And the last condition (50) is valid when gradients are small enough: 
IP. 



k 



N 



\P 



» 



a6 



N 



(51) 



So the condition (48) is always valid, condition (49) is valid when we investigate the 
instabilities near the perpendicular to the magnetic field and the condition (51 ) becomes 
only critical limitation for the approximation (47) of initial formula (46). 



4-.2. Ionospheric parameters 

To create a correct dispersion relation for heights 80-200km we should choose the correct 
approximations for ionospheric plasma. The most important plasma parameters are 
thermal velocities, hyrofrequencies and frequencies of collisions with neutrals. These 
approximate parameters are shown at the Table 1 (calculated for mid-latitude ionosphere 
using models MSIS, IGRF and IRI) 
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We also suggest that wavenumbers are within 0.1-lOm (sounding frequencies 15- 
1500MHz), and drift velocities \^ ao do not exceed 3000m/s. 

It is clear, that at altitudes 80-200km the following approximations are valid: 
\J 'm~ (vV^) << - electron hyrofrequency much higher than their thermal speed; 
("^^eoV^) << - electron hyrofrequency much higher than their average speed (even 
in disturbed conditions); ("^eo) << ~ electron thermal speed much higher than 
their average speed; v^:^ » electron- neutral collision frequency is much higher 
than ion-neutral one; We also use weakly inhomogeneous ionosphere approximation: 

Below there is a list of traditional approximations that are valid for E-layer 
[ni m [19], but is not vahd for the whole region 80-200km: (v^aoV^) << ^an - is 
not valid for ions at altitudes above 120km, not valid for electrons at heights above 
180km under very disturbed conditions; (l^ioVV^) << " is not valid for disturbed 
conditions; ili « - is not valid above 120km; >> v^^^ - is not valid at and 
below 80km. It also must be noted that during high disturbances the effective ion- 
neutral collision frequency v\^^ , that is used in dispersion relation, becomes dependent 



>> 



on electron density gradient and average ions velocity (38) and might be increased (or 
decreased) depending on the ion motion direction and gradients. It should be also 
noted that for high velocities the effective ion-neutral collision frequency v\^^ at heights 
approximately above 140-160km can become zero or negative. So, in this case the 
traditional approximation [H] of low Doppler drifts \u)c\ « WanNl is also invalid. 



4-3. Dispersion relation for weak gradients 

By substituting ( [39|43[44[ ) into ( [47) ) and after neglecting same non-zero multipliers: 



—T 



k Nb 



^1 + 



N 



(52) 



C2a {i^aN) = Z^eN, 



where 



)^ (Re(t 



N, 



(53) 



is hyrofrequency. 

Let us define new index f3 to describe another charged component: f3 = e,i; (3 ^ a. 
After defining the new parameters 

SoJ/SaN = ^PN - ^^aN = ("^aO - V' ^o) - il^ n) (54) 
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't = -Re{tN) = - (55) 

we will obtain the relations for second charged component as a function of the same 
parameter UaN- 



kh 



'^1 + 



(56) 



Below we will analyze the dispersion relation (29) in form: 
fix) = 

where 

fix = UaN) = ZpCiaix)C2pix) + ZaCi;3(x)C2a(x) 



(57) 



(5J 



By substituting ( 52p6 ) into (58), after neglecting non-zero multiplier, for single- 
charged ions iZi = —1) (most frequent approximation in this region of altitudes) the 
dispersion relation becomes the final one: 



fix) 



\X — IV, 

'2 ' t,^l 



IX 



Qfj (b't) + (v^f^^^ + iix + 5ujpaN, 



(t)> 



'l3nN 



+ iix + 6uj0a 



^ni{htY+{v'j:,^+ix)\-^Y) + 



(59) 



- iTa + 



bk 



fi^ (bVy + (^Z/J'^^ + iix + 5u)f3aN)y 



where 



X 



dip 



N y aO 



(60) 



and other parameters are defined by ( [35p8p4p5| and by solution of the zero-order 
approximation 

From obtained solution xq of dispersion relation (59) for given altitude dependence 
of the parameters one can always obtain the actual irregularity frequencies and 
decrements using relation: 

dijj 
~dt 



Xo 



(61) 
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From the dispersion relation 
approximation the presence of gradients H n of electron density logarithm change 
decrement (imaginary part of changes effective collision frequencies for ions 

with neutrals and makes them anisotropic at high altitudes. All these changes are 
proportional to the scalar product of the gradient of the electron density logarithm and 
average electron velocity. Actually this fact contradicts with current theories suggesting 
that in most cases only the electron density gradients perpendicular to the magnetic 
field must be taken into account [I9], since there could be conditions when ^ ao is not 
perpendicular to the magnetic field, for example in case of non-perpendicular magnetic 
and electric fields. 



59p0pT ) it becomes clear that in the first 



Low-frequency solution 



4-.4-1- The solution nearest to zero The dispersion relation (59) has 6 solutions, and 
in basic case all the solutions can be found only numerically. Lets find the simplest 
approximate solution - nearest to zero. The solution nearest to zero has a clear physical 
sence: in absence of average plasma drifts and gradients plasma can be supposed as 
static and irregularities should be static, i.e. 

In presence of weak drifts and gradients we can suppose that the solution is close 
to zero. 

To find the solution nearest to zero we will use zero order Newton solution (see, for 
example [21]): 



Xo = - 



\dx) 



(63) 



x=0 



By substituting the basic relations: 

/(O) = mp {6ujpo,n) [^^paN - iJ^f 



bk 



^1 + 



+ 



(64) 



"/3 



bk 
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+ i6uJpaN) 



+ 



n}(bk) +( 

+2imp {6u)paN) [^^PaN - it^finN) ' 



-2z (T„ + Tf, 



+ 



bk 



- iduJf^aN 



(65) 



The low-frequency branch ( [63|64|65[ ) is pretty complex so we will investigate it at 
different ionospheric heights. 



4-4-^- Obtaining traditional solution at 80-120km heights Let us analyze the branch 
p3|64|65| ) for the typical ionospheric heights 80-120km. Within standard for E-layer 



assumptions of x, SuieN « i^lnNj ^\nN t magnetized electrons and unmagnetized ions, 
and neglecting 5u)ieN « ^InNi obtain following equations for function and its first 
differential (neglecting in first differential by all the terms, proportional to SuJieN or 
Te + Tj, based on suggestion that Doppler shifts for ionacoustic or average velocities are 
sufficiently small in comparison with z^i^jv): 



m 



2 / + ,, \ 2 



t.fi 

inN 



{n + T,)(ki 



(66) 



dfjx) 
dx 



x=0 



-^m, (ut^^) {ut^^y (nl {bty + {ut%y (t) 
From ( 67||66 ) and Newton method (63) we obtain the solution, nearest to zero: 



(67) 



Xo 



1 + 







(''enjv)^ 


( t,p. \ t,^ / 


bk 


2 




(^enJv)M 



(6J 
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Considering (61) we obtain the following solution for plasma irregularities: 



(hp 
dt 



'I'+l 



— I- 



(69) 



I'+l 



/,tJ. t,fj. 



bk 



( V 



(^) (-i^i 



enN J 



(70) 



In gradient-free case the solution (69) looks exactly as the standard one [19], in 
presence of gradients the solution differs from the standard one, most probably due to 
custom direction of gradients and custom orientation of velocities. 



4.4-3. Fully magnetized case or instabilities at altitudes above 140km Let us analyze 
the branch ( [63p4p5| in case of sufficiently high altitudes, when both types of charged 
particles are magnetized (i.e. from about 130-140 km). At high altitudes we can neglect 
the difference in electron and ion velocities in comparison with their absolute values 
(both components are magnetized and move with almost the same velocities): 



\Re{x)\ » \6uJi3^ 



aN\ 



(71) 



Supposing dujpaN = 0, and when investigating wavevectors perpendicular to the 
magnetic field, the solution becomes simplier: 



/(0) = -(T„ + T; 



\yanN ) y^PnN ) 



+ 



-2i (r„ + Tp 



(72) 



(73) 



^l3nN^anN 



0nN + ^aniv) 



After simple arithmetic and by taking into account magnetized plasma and typical 
ionospheric conditions (a- electrons): 



^anN ) 'i^^P 



Ql » 



we have: 



/(0) = -(T, + T^ 



y^anN ) \'^PnN 



« V. 



anN 



ID 



IV, 



f3nN 



-2% (T„ + Tp 



( ■ t,fj. 



t,fi y 

2 



(74) 



(75) 



'^(3nN y^anN ) 
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Taking into account the typical ionospheric conditions: 



k 



''enN 



^enN 



Xo - 

dijj 
~dt 



t,fj, 

''enN 



m me 

{Te + T,] 



- i'^N^eO + ^eO ^ 



From (78) the condition for the growing solution Im < becomes: 



l^N'^eO > D 



where 



D 



V. 



Al. 



t,fj, (Te+T,) 



enN~ 



A± 



nii 



(76) 
(77) 

(78) 

(79) 
(80) 



is the so called coefficient of ambipolar diffusion [T5] . 

For typical ionospheric conditions the growth condition ( 79 ) can be estimated as: 

2 , , 



k 



The spectral offset for these irregularities (78) is exactly the Doppler drift in crossed 
fields (and defined by zero-order solution (|4]|5|): 

' dip\ 



Re 



dt y 



2) 



It is necessary to note that the solution is obtained in weak gradients approximation 



(51 ) valid when: 



k 



^enN 



» 



N 



(83) 



It should be noted that possible relation of the ambipolar diffusion with 
irregularities existence at these heights has been noted in [18j, but the problem was 
not investigated in detail. The close condition for irregularities growth at high altitudes 
Ve > TjDAi. was also obtained by (TU], but with another proportionality coefficient r] and 
in qualitative analysis of a model case. 

Due to the solution (78) corresponds to the same branch ( [63|64|65[ ) as well-known 
gradint-drift instabilities ( 69||70 ) , below we will call this kind of solition as 'fully 
magnetized gradient- drift instabilities' (FMGD) to stress that this is the same gradient- 
drift branch but in a bit different conditions. 

Starting from early 1960s [3l|2] at equatorial HF and UHF radars researchers observe 
an unique type of echo, the so called 150km equatorial one. There are some theories 
to explain it (for example [20l [221 ESI [TBI |9]) , but the exact physical mechanism of 
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it is still unclear [El El [7j. The geometry at equator (horizontal magnetic field, almost 
upward drift velocity) allows us to use standard upward vertical gradient as a source 
for generation of this kind of instabilities. For sounding frequency of equatorial radar 
Jicamarca {k ~ 0.3) and for ~ 10~^ 10~^[m~^] (standard vertical electron density 



gradients) the growth condition (79) becomes: 



l^eo > 20-^200 [m/s] (84) 



It should be noted that for gradients higher than (83) {K^ > 10 ^) one should take 



into account all the terms in (46) instead of using only (47), but this should not affect 
too much the observed effect. 

To analyze the properties of the echo, some modelling has been done using the latest 
Internation Refference Ionosphere (IRI-2007) model. The height and time dependence 
of the electron density gradients are most important for the generation of this type of 
irregularities. We have analyzed 13 years period (1990-2002) using IRI model (for typical 
non-disturbed conditions /10.7 = 150, Ap = 10) and obtained the following results. 

At Fig.l the altitudinal dependence of K]sr is shown. Points are the hourly values 
over the whole period of 13 years. As one can see, there is a maximum Kn > 10^^ at 
heights 135-180km. So, this kind of irregularities could arise at heights 135-180km, and 
it corresponds well with the experimental observations statistics O |8] . 

At Fig. 2 an hourly dependence of the Kn is shown, as a function of UT for heights 
140-200km. As one can see, the time dependence of the gradients has a most intensive 
maximum between 14:00 and 19:00 UT (9:00-14:00LT). This also corresponds well with 
the experimental observations 



The dependence of irregularities frequency (82) corresponds well with the empirical 
models [6l [21] and allows to interpret the experimental data as Doppler frequency offset 
due to electron drift in crossed fields. 

According to the experimental observations, the echo starts with Ve > 10m/ s, 
according to our calculations it should start with Vg > 20 ^ 200m/s. One of the 
mechanisms allowing to lower the speed hmit was suggested in [20]. They suggest that 
acoustic-gravitational waves can be responsible for the triggering the instabilities. In our 
terms, the acoustic-gravitational waves will produce gradients more than > 10~^, 
and this will produce this type of irregularities even at lower velocities, for example at 
Ve > 10m/ s. Another possible mechanism that will lower the velocities necessary for 
generation of the instability is an observation of high step-like gradients at these heights 
from the rocket data (see for example [21] )• They should also produce the increase of 
Kn high enough for lowering the speed limit. 

Summarizing all said above we can suggest that the FMGD instabilities can be 
the source of 150km equatorial echo and this theory can be used for experiment 
interpretation. 
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Table 1. Thermal velosities, hyrofrcquencies and frequencies of collisions with neutrals 
for typical ionospheric conditions. Frequencies in Hz, Velocities in m/sec, height in km 
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5. Conclusions 

In the paper within the approximation of the two-fluid magnetohydrodynamics and 



geometrooptical approximation the dispersion relation (59 60 35 38) at 80-200km 



altitudes was obtained. The relation describes ionacoustic instabilities of the ionospheric 
plasma at 80-200km altitudes in three-dimensional weakly irregular ionosphere. 

It was shown that not only electron density gradients perpendicular to the magnetic 
fleld should be taken into account when investiagting ionospheric instabilities, but 



gradients along the average drift velocity (59, 60, 36, 38). 



The dispersion relation obtained has a form of the 6-th order polynomial for the 
oscillation frequency. 

It is shown, that a solution branch exists that grows with time and describe 
instabilities both at 80-120km heights and 135-180km heights. 

For altitudes 80-120km the solution close to the standard one (69 70) and 
corresponds to the Farley-Buneman and gradient- drift instabilties. 

The difference between obtained ( [63|64[65[ ) and standard solutions [19] becomes 
signiflcant at altitudes above 140 km, where standard one is not valid. As the analysis 



shown at these altitudes the solution grows with time (78, 79). The conditions for the 



growth is the presence of co-directed electron density gradients and electron drifts and 
perpendicularity of line-of-sight to the magnetic fleld. These conditions are regularly 



satisfled at magnetic equator for expected conditions (84). Detailed analysis has shown 



that this solution could explain a lot of properties of 150 km equatorial radioecho - the 
ionospheric phenomena that has no explanation for more than 40 years. 
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Figure 1. Km dependence on height over the 1990-2002 years, points are the hourly 
values. 



Appendix A. Approximations used 

The theory hmitations are hsted, mostly following to the [T71 |15l [H] . 

At the altitudes 80-200km we suggest that the following conditions are satisfied: 
« VLi - all the basic plasma parameters has only slow variations and plasma 
supposed to be quasistatic; 

Sen << 1 - average loss of energy of electrons with neutrals is small enough; 

z/jj << z/j„ - ion- ion collisions are rare enough to take into account only ion-neutral 
collisions. Not valid above 200km. 

i^ei ~ ^'ee << ^en " elcctrou-ion and electron-electron collisions are rare enough to 
take into account only electron-neutral collisions. Not valid above 200km. 

k±Pe << 1 electron hyroradius much smaller than wavelength. 

Pe « Ae << k^^ « L\\ - plasma is quasihomogeneous enough for GO 
approximation to be valid. 

Xdk << 1 - wavelength is much bigger than Debye radius. 

Seii^ei « i^in - ueccssary for independent thermalization of ions and electrons, in 
this approximation the average collision frequency does not depend on particles velocity 
or motion direction jlTj. 

Von = - average speed of neutrals is much smaller than electrons and ions speed. 
Allows us to neglect neutral motions. 
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Figure 2. K^q dependence on time at 140-200km over the 1990-2002 years, points are 
the hourly values. 

In more details the MHD validity conditions can be found in [T71 [15] . 
Appendix B. Inversion of the matrix operator Pq, q (IAQV approximation) 



Lets analyze inversion of the matrix operator P^ q (22 25). One can see, that in special 
case P7 = the inversion is very easy. In this case by taking into account that 'P q is 
static, we can create the coordinate system, based on unity vector h = Pq/\Pq\, which 
is antiparallel to the magnetic field. In this case we can write: 



po7 = a7ii + p^y^ + pji. X h 

Where 1 1 , _L means parallel and perpendicular to the h. 



(B.l) 



By making scalar and vector products of (B.l) with h we have: 

h P^h X (7^ X h) 



h{P,f) = P,hf\\ 

h X {Pj) = P,(b X 7 



(B.2) 



From first equation (|B.2|): 

7,1 



MPo f 



p. 



(B.3) 
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After comparing the second equation in (B.2) and its vector product with b we 



have: 



b X (b x'f 



&x(Po/)-P5(feX /^) 
bxibxjPo f ))-P6bx{bx{ f ^Xb)) ^l^^^ 



Therefore, by taking into account the properties of double vector product: 

bx(Po7)-P5(fcx7^) ^ 



bx(fex(Po / ))-P6(feX /^) 



So 



(bx f 



nfex (Po/)+P6(&X (&X (Pq/))) 

Pi + Pi 



and (after making vector product with b and some vector algebra): 
^ Pejb X (Po7)) - P56 x(bx (Po7)) 



Summarizing (B.3) and (B.7) we have / = f ± + f 



_ P,P,{b X (Po7)) - Pjb xjbx (Po7)) 

^ (^1 + Pi) A 



+ 



+ 



(p| + p|)6(6(Po7: 



(^1 



Pi) A 



or 



b(bf)Pi + Pif +P,P,{bx f 



P,{Rl + Rl) 
It is clear that this approximation is valid when: 



Pa57 + 7 X 



>> 



P^rJl 



Qualitatively one can estimate the orders of terms: 



^la >> VaoLy^ 



where 



ay 



(B.4) 



(B.5) 



(B.6) 



(B.7) 



(B.8) 



(B.9) 

(B.IO) 
(B.ll) 

(B.12) 



For maximal ionospheric disturbances up to 200km height we can estimate Vao < 
3000[m/s], > lOOHz, > lOHz, > lOOHz. In this case the validity condition 
has the form 



Lv» 300 [m] 



(B.13) 
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Summarizing, in very disturbed ionosphere the characteristic changes of parameters 
should not exceed couple hundreds meters, for less disturbed conditions these limitations 
becomes even weaker. So the obtained approximation for P~J^ (B.9) is valid for most 
part of cases below 200km. 
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